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Abstract 

The S wave pion deuteron scattering length is presented in an effective field theory 
approach of the two-nucleon system. We include pions as dynamical particles which al- 
lows us to calculate pion re-scattering contributions inside the deuteron. Two-nucleon 
two-pion contact interactions with unknown parameters not determined by chiral sym- 
metry need to be introduced in order to renormalise the appearing divergences. By 
choosing their values accordingly we are able to accommodate the available experi- 
mental data. 
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1 Motivation 



A systematic framework to calculate low-energy scattering processes of hadrons is provided 
by Chiral Perturbation Theory (^PT), the effective field theory of QCD at low energies. For 
example, he determination of the pion nucleon scattering length has drawn a lot of attention 
and still remains a topic of interest. The chiral corrections to the Weinberg- To mozawa 
current algebra theorem according to which the iso-scalar (iso-spin even) S wave pion nucleon 
scattering length := \{aT^p + flTm) is zero at leading order (LO) have been calculated by 
Bernard, Kaiser and MeiBner 0|. Their next-to-leading order (NLO) prediction has 
recently been refined by determining the unknown parameters in a comparison of non- 
zero energy vrA^ scattering at the chiral orders and to Koch's Karlsruhe 0, 
Matsinos' EM98 and the VPI/GW group's SP98 0] partial wave analyses. The obtained 
range OhbxPT — [~0.01 . . . -|- 0.006] is still compatible with zero. Recently, Gasser et 
al. pointed out that electro-magnetic iso-spin breaking effects could provide corrections of 
up to 7% to the LO result. Because a+ is very small, such effects are stronger than the naive 
expectation that they are naturally suppressed by powers of the fine structure constant a = 
1/137, contributing at most 1%. Another recent HB^PT analysis of pion nucleon scattering 
to order attempting to connect threshold and near-threshold parameters with the low 
energy theorems of chiral symmetry found on the basis of [|5| that = —0.008 0. 
Direct extractions of a"*" from the same ttN partial wave analyses yield, on the other hand, 
(+0.0041 ± 0.0009) g and +0.002 §. 

On the experimental side, the best extraction of a+ comes from measurements of elastic 
scattering and single charge exchange in pionic atoms. In contradistinction to phase shift 
analyses, no extrapolation to zero energy scattering is necessary, and electro-magnetic effects 



are considered small. Most recently, the ETHZ-Neuchatel-PSI [jlO| collaboration found 
a^p = (—0.0022 ± 0.0043) directly from the line shift and width change in pionic 
hydrogen, assuming iso-spin symmetry. This is compatible with zero and with the HB^PT 
result, but with smaller error bars. The ongoing experiment R-98.01 at the Paul Scherer 



Institute aims to reduce the error by an order of magnitude |11 
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The value of the pion deuteron scattering length, which was measured by the same 
method as for pionic hydrogen to be a^d = [(-0.0261 ±0.0005) + i(0.0063±0.0007)] m^i p 



has also been used to constrain . However, since - as mentioned above - the deuteron is 
not a purely iso-scalar nucleon target, binding and especially pion re-scattering effects be- 
tween the two nucleons have to be accounted for. Conventional potential model approaches 
which combine deuteron and hydrogen data have been utilised: Baru and Kudryatsev, e.g., 
used multiple scattering methods and quote a value o^^l^^ bk ~ (— 0.0015±0.0009) [1131 , 
however, their errors are substantially underestimated as shown in [^. An analysis by 
Landau and Thomas [jT3| gives a^-^^^ lt ~ (+0.0016 ± 0.0013) which recently has been 
updated by Ericson, Loiseau and Thomas [Q to flphcn, elt — 

(-0.0017 ± 0.0002(stat) ± 

0.0008(sys)) Tn~^ . In that work it was also mentioned that is also a major source of un- 



certainty in the extraction of the coupling constant Qt^nn [0]- With the negative value for 
a"*" given in [jl^, S'^attv/I^tt) = 14.17±0.20 is extracted in contrast to most other approaches 



which favour qInnH.^'^) ^ 13.7 ± 0.1, see Table lin [0. 

Given the unsatisfactory experimental situation for it seems preferable to take into 
account binding and pion re-scattering effects between two nucleons in a more model- 
independent fashion, in order to constrain from pion deuteron scattering. Such a frame- 
work is provided by xPT which incorporates the chiral symmetry of low-energy QCD. Em- 
ploying the so-called Weinberg counting [|l^ in which the pions are treated non-perturbatively 
and iterated to infinite order in ladder exchange diagrams, Beane et al. calculated 
a-j^d- Their result for a-j^d does not involve any undetermined parameters, however, the non- 
perturbative effects responsible for nuclear binding are accounted for using phenomenological 
deuteron wavefunctions which introduces an inevitable model dependence. It is not clear 
to what extent the phenomenologically based deuteron wavefunctions employed in [T7| are 
constrained from chiral symmetry, and one may therefore pose the question: 
Does chiral symmetry predict the pion deuteron scattering length or is additional phe- 
nomenological input needed in order to fix new unknown parameters? 
Clearly, this cannot be answered within the Weinberg framework which always employs 
deuteron wavefunctions in order to account for the non-perturbative behaviour of pion ex- 
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change. One must resort to a different approach which is given by the method of Kaplan, 
Savage and Wise [|1^ (ENT(KSW)) where, in contradistinction to the Weinberg scheme, 
pion exchange between the two nucleons is treated perturbatively. It has been illustrated in 
a number of papers that ENT(KSW) calculations at next-to-leading order (NLO) are - for 
small momenta - in good agreement with experimental data. At higher orders, however, the 
perturbative expansion is expected to fail in the deuteron channel |T^. Since our primary 
interest is to clarify the existence and importance of contributions to the pion deuteron scat- 
tering length not fixed by chiral symmetry, we will restrict ourselves to a next-to-leading 
calculation in the KSW scheme. As we will see, this issue can already be investigated at 
next-to-leading order so that we are not concerned with problems arising at higher orders. 

The article is organised as follows: We start by presenting the Lagrangean (Sect. 
In Sect. ^, the vrrf scattering amplitude at zero momentum and its scattering length is 
presented, followed by a discussion of our findings in Sect. including the comparison to 



the EFT calculation of Beane et al. [1? . The final Section contains our conclusions. Two 
Appendices summarise details of the calculation and renormalisation procedure. 



2 The Lagrangean 



We now present the terms of the most general chirally invariant Lagrangean consisting 
of contact interactions between non-relativistic nucleons, and between nucleons and pions, 
which are relevant for our NLO calculation. 

The pertinent terms satisfying the QCD symmetries are in the zero and one nucleon 
sector (see e.g. pO[] ) 



C.,.N = ^tr[(D^St)(Z}^S)]+ ^ cc.tr [7Wq(St + S)] + 
o 4 



2M 



)N + QAN^A-aN + 



(2.1) 



2u; ci tr[A^q(S + S^)] + A { C2 - ^ ) Al + Ac^ A^A^ 



8M 



N + ... 



where = (^) is the nucleon doublet of two-component spinors, M = 938.918 MeV the 
iso-scalar nucleon mass, and a (r) the Pauli matrices acting in spin (iso-spin) space. The 
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field ^ describes the relativistic pion, for wliicli we clioose tlie sigma gauge for convenience, 



pa 



= e{x) = ^1 - 2 + iv^ . (2.2) 

Dfj^ is the chirally covariant derivative = + Vf^, and the vector and axial currents are 

V, = l{^d4^ + ed,0 , A^='-{^d,e-ed,0 ■ (2.3) 

The pion decay constant is normalised to be /tt = 130 MeV, qa = 1.27, Aiq = diag(m, rh) is 
the quark mass matrix where we work in the iso-spin limit m„ = rrid = rh, and the constant 
u is chosen such that = 2ujrh at the order examined here, where we use the iso-scalar 
value m^r = 138.039 MeV for the pion mass. The coefficients ci, C2, C3 encode high-energy 
physics integrated out in HBxPT and need at present to be determined by experiment. 
Since at the scales considered, the momenta of the nucleons are small compared to their rest 
mass, the nucleons are treated non-relativistically at leading order in the velocity expansion, 
with relativistic corrections systematically included at higher orders. Thus, the relativistic 
HBxPT Lagrangean is reduced to the form shown above. 

The germane terms in the two nucleon Lagrangean are (see also ||19|| ) 

Cnn,.nn = -Co (iVTpW)t (iV^PW) + ^ \{N^P'Ny (N^P^D - DfN) + H.c. 

o L 

- I A tr[A^q(S + T))] {N^P'Ny (N^P'N) + (2.4) 
+ (^E2 ti[Al] - F2 tr[l2] j (N^p'N)^ (N^P'N) + ... , 

where is the projector onto the S wave of the iso-scalar-vector channel, 

= {r,)i , (2.5) 



and the parameters of the Lagrangian Ci,D2, E2 and F2 are not constrained by chiral sym- 
metry, but can be extracted from experimental input or estimated with additional model 
dependent assumptions. E.g., Cq can be related to the binding energy of the deuteron and 
C2 to the nucleon-nucleon scattering length [|l^. The remaining parameters D2, E2 and F2, 
on the other hand, have not been determined yet. 
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3 Calculation 

The pion deuteron scattering length a^^ follows from the amplitude A-^d at zero momentum, 

1 r Tn 1 

which in turn is decomposed into a contribution in which the pion scatters off only one 
nucleon (Fig. 0), and a term with two nucleon interactions (Fig. 

A. = + A'^f'^^ . (3.2) 

The first amplitude starts at LO, 0{Q^) after wave function renormalisation, while the 
second one is NLO, 0{Q^). 

3.1 One Body Contributions 

We first discuss briefly how to embed the well known HB^PT result for the iso-scalar S wave 
pion nucleon scattering amplitude A~^ into the deuteron. The one body contributions to nd 
scattering as shown in Fig. [l| consist of the LO and NLO iso-scalar pion nucleon amplitude 
A'^ at zero momentum (diagrams (Qa) and (|1^)), corrections arising from wave function 
renormalisation and NLO deuteron effects (diagrams (0c)). Since the latter diagrams can 
be absorbed into a re-definition of the deuteron source used, it is no surprise that their 
contribution to A-j^d at zero momentum cancels with the NLO wave function renormalisation 
of the LO amplitude. 

As the pion can scatter off the proton or neutron inside the deuteron, the one body 
amplitude is twice the physical amplitude A'^ only, independent of any deuteron observable. 
This is formally confirmed in an exemplary calculation in App. ^ Relativistic effects do 
not enter at the order we calculate. Therefore, 



A Ibody 



Am, 



2 



J 77 



2A+ (3.3) 



2 8M QA-Kfl 

with the bare quantity Abare(/i) := 2 (c2,bare(/i) + C3,barc(/u) - 2ci_barc(yu))- Notice that both 
graphs {hi) and (62) in Fig. |l] depend on the renormalisation scale as the nested integral 
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(m^,0) 



A+ 



Cl,C2,C3 



i-Bfi) 




(-5,0) 




+ 



(a) 



(61) 



(62) 



+ 




(cl) 



Figure 1: The contributions from : (a); LO; {b): NLO in A'^ ; (c); NLO deuteron cor- 
rections. Shaded circles denote the deuteron source. Graphs obtained by permuting vertices 
or external lines are not displayed. 

is linearly divergent. In the MS or MS scheme usually chosen in HBxPT, the divergence is 
discarded (Abare(/^) = A'^^), but it is manifest in the PDS scheme which for consistency has 
of course to be used also in the one nucleon part of the calculation. Therefore, the HBxPT 
parameters Ci, C2, C3 start to depend on the renormalisation procedure at NLO and are not 
observables. It is however obvious that the regulator dependence is absorbed into the NLO, 
0{Q) part of the combination of coefficients A(//) = A*^'^^ + A^^\fi), in agreement with the 
power counting. The LO part, A'^^-', is fi independent. The value for A in the MS scheme 

2 

is hence easily translated into the PDS scheme as Abare(/^) = A'^^ + ^|^j2 , making the full, 
physical one body amplitude (|3.3| ) exphcitly fj, independent. 

Recently, updated values for the parameters Cj have been determined in an 0{Q^) and 
HBxPT fit to three finite energy pion nucleon scattering analyses 0, §]. As already 
noted in Refs. the iso-scalar S wave scattering length 



a 



+ 




A+ 



(3.4) 



can however not be determined precisely because of a numerical cancellation which may 
signal physics at small scales not yet understood. The HBxPT analysis predicts the range P, 
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i 

a+3^p^ = [-0.01 . . . + 0.006] m-^ = (-0.002 ± 0.008) m^^ , (3.5) 

compatible with the various experimental and phenomenological extractions discussed in the 
Introduction. The comparatively large range comes from the use of different partial wave 
analyses for the vrA^ amplitudes; the theoretical uncertainty from HBxPT is considerably 
smaller. 



3.2 Two Body Contributions 

Processes involving deuteron correlations, Fig. enter at NLO, 0{Q^). The sum of the 
two diagrams Fig. |] (al) and (a2) is independent of the parametrisation of the pion field. 
The two pion one nucleon vertices in (b) stem from the chirally covariant form of the time 
derivative in the kinetic energy term for the nucleon in (|2.1|). 




Figure 2: The two body contributions from pion re-scattering [a, b) and two body contact 
interactions (c) at NLO. Graphs obtained by permuting vertices or external lines are not 
displayed. 



The diagrams with a pion in the intermediate state, Figs. H (a) and (6), contain log- 
arithmic divergences. In order to see how they manifest themselves in a more traditional 
setting, consider the graphs when a deuteron wave function is used instead of the point-like 
source with the correct quantum numbers which is the starting point of the standard field 
theoretical treatment. In this case, the contributions are proportional to the expectation 
values of two operators [|r^ ^ 



A 



2body{a) ■ ■ ^2) . 



9 



^T^l^ oc (1) (3.6) 

between deuteron wave functions, where q is the momentum transfer between the nucleons. 
Thus, the graphs seem to probe the deuteron only at momenta of the order of g ~ m^r 
in Fig. 1^ (a) and g ~ in Fig. ^ ih). The Fourier transformation into position space 
reveals however that the expectation values of the operators ^ and ° are probed, testing 
the deuteron wave function at arbitrarily short distance. Usually, the UV part of the 
deuteron wave function is parameterised by adding to the long range pionic potentials strong 
phenomenological terms mimicking short range repulsion. The deuteron wave function 
is thus suppressed at large momentum transfer at the price of adding some dependence 
on the shape and parameters of the short distance potential. Thus, the result for these 
amplitudes depends also on short distance physics subsumed into an arbitrary, unphysical 
cut-off parameter. 

It is therefore no surprise that a logarithmic divergence appears in the pion exchange 



diagrams. The necessary integrals are tabulated in ||2^, so that the bare amplitudes are 

+ 27 



.2body(a) _ g\'ml-i _ ...^ _ 

•A.,, bare " ^^^f, ^^^37 

.2body(6) _ _ /^r - 91n 

■^■Kd, bare " ^ J4 ' ^ ^ 



(3.7) 

/i 

with a divergence F := F[4 — d\ + + 1 in c? ^ 4 space-time dimensions Q Therefore, 
the bare two pion two nucleon contact interactions Fig. |^ (cl) and (c2) entering at the same 
order are also necessary in order to consistently remove all regulator dependence in the 
total, physical two body amplitude. 

■4^:1? = ^ [D2, barc(/i) + E,^ barc(/i)] (7 ' f^)' (3-8) 

They serve the same purpose as suppressing the deuteron wave function at short distance by 
a cut-off. As must be expected from the power counting, the combinations -D2, bare(/^) {l^l^Y 
and £'2, bare(/^) (7 — scale as . The strengths of both two pion two nucleon contact 
^Following [2^ , we prefer to work with a definition of F whieh is independent of the scale /i. in contradis- 
tinction to e.g. [pO[. 
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terms cannot be predicted at the present time. Albeit the parameter D2 is encountered as 
the chiral symmetry breaking contribution to NN scattering at NLO, in this process its 
renormahsed value cannot be dis-entangled from the parameter Cq^^ which respects chiral 



symmetry [T^. The strength E2 is not fixed, either. Unfortunately, we are also not able to 
extract their values from additional experimental input so that we cannot predict the pion 
deuteron scattering length. 

Nonetheless, EFT allows to estimate the natural size of the combination D2 + E2 after 
specifying a prescription to handle the divergences in three and four dimensions. Only after 
renormalisation is one justified to compare the sizes of the various contributions. In the 
PDS scheme, the theory is usually renormahsed (and the power counting made manifest 
even before renormalisation) by choosing the renormalisation scale to be natural |]I8[: 



/i — > , r — i> 2 (3.9) 

With the latter choice, the divergences in 4 dimensions are removed by demanding that 
the pionic contribution to zero momentum scattering between two nucleons disappears, see 
([B.lp in App. |B|. Power counting and dimensional analysis dictate then that 

D2 + E2 = -^-^-^, (3.10) 
M Ann (At- 7) 

where the magnitude of the dimensionless parameter is of order unity, \Z\ ^ 1, if the 
naturalness assumption holds. Its sign remains undetermined. The scale Ann ^ 300 MeV 
enters in the chiral expansion for NN scattering |TB| and subsumes al short distance physics, 
i.e. all effects of particles not contained as explicit degrees of freedom in the Lagrangean, 
like the A, p meson exchange etc. 

However, it must be stressed that in the approach taken here, a decomposition of the 
physically observable nd scattering length into parts related to the pion re-scattering di- 
agrams of Fig. (a) and (b) separately is strictly speaking meaningless: None of these 
diagrams is renormalisation group invariant (i.e. cut-off independent) on its own, and only 
combinations of these diagrams with the two pion two nucleon contact diagrams of Fig. |^ 
(c) form observables free of divergences, i.e. independent of fi and F. After renormalisation 
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performed in App. 0, the scattering amplitude reads 

1 



2 

^2body _ "^7r7 



2 In 


[27] 


9\l 




A* . 


6 I 



+ In 



27)' 



A* 



(3.11) 



and contains only one un-determined, physical parameter A*. This dimension-ful number 
parametrises the renormalisation group invariant strength of the contact interactions D2 
and E2 between nucleons and subsumes effects from the deuteron wave function at short 
distances. It needs to be determined from experiment or from a microscopic calculation of 
A^A^ scattering in QCD. A* is expected to be of the order of the natural low energy scale 
(m^ or Aqcd) since all dependence on higher scales integrated out has disappeared with 
renormalisation. 

As mentioned above, relativistic effects do not enter at the order we are working. Nev- 
ertheless, it may be worthwhile to comment on their contributions. Relativistic corrections 
to the energy-momentum relation are accounted for by inserting perturbatively higher di- 
mension operators, the lowest one being p^/(8M^). They start at N2L0 and are suppressed 
by factors of the nucleon mass, and not by K^n ~ 300 MeV. Such effects are therefore 
small compared to other corrections which enter formally at the same order N2L0. Our 
calculation only includes NLO effects. For our purposes, they are therefore negligible, and 
may be regarded to lie within the given error bars. 



4 Results and Discussion 

With the data at hand, we cannot predict either the iso-scalar or the pion deuteron scattering 
length in a unique way from our calculation due to the unknown physical scale A* stemming 
from the combination D2 + E2 of counter terms. These unknown strengths of the two-pion 
two-nucleon couplings enter already at NLO, i.e. at the same order as pion re-scattering. 
In order to determine their sizes, we use as inputs the experimental value for the pion 
deuteron scattering length and a value for . Clearly, requiring self-consistency forbids to 
consider phenomenological extractions of a+ into which Qt^cL entered. We use as choice either 
the value from the pionic hydrogen experiment, = (—0.0022 ± 0.0043) |T^, or 
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the HB^PT prediction, Eq. ( p.5|) 0, Both are compatible with zero, and the HB^PT 
result has an error bar accommodating also the phenomenological partial wave analysis 
extractions ^, 0| cited in the Introduction. We find 

A* = (26211^^) MeV when fitted to a+p^^^ 

A* = (269+^70) MeV when fitted to Oj^bxPT § > (4-1) 

which contains a large error bar from the uncertainty in a"*". The value of A* increases as 
a"*" increases. The Logarithms in ( |3.11| ) are indeed of order 1. 



We now consider the sizes of two renormalisation group invariant subsets constructed 



out of the two body scattering result (|3.11|) : The first one is the term independent of 



qa- The second one is quadratic in qa, and is easily estimated to be suppressed by a 
factor 1/12. We reproduce the well known result that the dominant pion re-scattering 
contribution stems from physics unchanged by taking the chiral limit, and A* ^ A^^. Albeit 
the definition of A* depends also on qa ( |B.8| ), this polynomial separation with respect to 
Qa is at present a good analogue to the values quoted in the Weinberg approach for the 
pion re-scattering contributions. Fig. ^ (a) and (b). Although much simpler, our results for 
these renormalisation group invariant combinations are close to the findings of the Weinberg 
approach, see Table [l|. This should not be too much a surprise since it was shown above 
that the pion re-scattering diagrams test after renormalisation only the tail of the deuteron 
wave function, at momenta not larger than m^. Figure ^ depicts the dependence of the two 
body contributions to the Tid scattering length on A*. We finally summarise all one and two 
body contributions to ajrd in Table || for the extraction using a^p^^. 

A similar analysis can also be made for D2 + E2. Taking a^J^p from the pionic hydrogen 
experiment, we find using the prescription ( |3.10| ) 



2 = 1.0^1.5 (4.2) 

with a large uncertainty due to the error bar in the value of a^p For the central value Z = 
1.0 the contributions from the individual diagrams are a^f^^^""^ = +0.0005 m~^, a^f'^^^'''^ = 
—0.028 ?Ti~^, a^J^'^^^^^ = +0.006 We note that choosing a different prescription for 
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approach 






KSW , RG invariant 


-0.0004 


-0.0210 


Weinberg 


-0.0007 ±0.0002 


-0.0196 ±0.0005 



Table 1: Comparison of two body corrections to a-j^d in- the KSW power counting with A* 
determined from and in the hybrid approach based on Weinberg's power counting 

taken from Error bars from the uncertainty in omitted. The error bars in the 

Weinberg approach stem from using different model wave functions for the deuteron. 
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Figure 3: Sizes of the renormalisation group invariant contributions to aT^d, ( \3.11\ ), depending 
on the physical parameter A* . Dashed line: contribution binomial in qa; dash-dotted line: 
contribution independent of qa! solid line: total. 

removing the 4 dimensional UV divergences, F — > 1, leads to pion re-scattering contributions 
numerically deviating by less than 8% from those obtained by Beane et al. |1^: a^f^^ (F 
1) = -0.0008 m-\ af/^^'^(F 1) = -0.019 m-\ af/^^'^(F 1) = -0.002 with 
Z{r 1) = —0.4 =1= 1.5. As the error bars indicate, the magnitude of D2 + E2 increases 



as a"*" increases. Again, the contribution from the counter terms is clearly not unnaturally 



large and is roughly 10 — 20% of the dominating pion re-scattering amplitude a^^ 



2body (b) 



This 



implies that a 10% uncertainty in the parameters D2 and E2 yields an error of about 1% 
in the total scattering length. The effect of the counter terms is thus even smaller than the 
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diagram 


contribution to a^^^ ["^tt^] 


1 body from a+p,^^ 


-0.005 ±0.009 


2 body, total 

2 body, g\ 
2 body, ^2 


-0.021^0.009 

-0.021 T 0.009 
-0.0004 ±0.0014 


exp 


-0.0261 ±0.0005 



Table 2: One and two body corrections to 0,^0^, choosing as input the experimental extraction 
of from pionic hydrogen experiments / [7^/ . 

power counting suggests. 

One may also assume that the two body counter terms D2 and E2 are saturated by a 
mechanism which the phenomenological extraction of a+ from the pion deuteron scattering 
length by Ericson et al. [|14| can capture correctly. Taking their value a^^^^^ ^^rj., we obtain 
A* ^ 280 MeV, Z ^ 0.8. However, such an approach violates the spirit of self-consistency 
at the basis of our calculation. 

5 Conclusions 

We presented a calculation of the pion deuteron scattering length in Effective Field Theory 
with perturbative pions. In this scheme proposed by Kaplan, Savage and Wise, knowledge 
of the iso-scalar pion nucleon scattering length a"*" does not suffice to determine ajrd directly 
due to two unknown short distance effective interactions coupling two nucleons to two pions. 
Their strengths D2 and E2 can be subsumed into one physical unknown which we estimated 
in the range A* ~ 270 MeV, but with sizeable error bars. Employing A* ~ 270 MeV, we 
can reproduce the numerical results for the pion re-scattering graphs in the work by Beane 
et al. |l^ who use the Weinberg scheme. Unfortunately, the experimental error on is 
not small enough to constrain A* (or D2 + E2) to a more precise value. The 10% accuracy 
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assigned to our NLO calculation is not only comparable to the uncertainty induced by 
experiment in the sizes of our counter terms. It also masks iso-spin breaking contributions 
to a"*", which were shown in ^ not to exceed 7% of the pion nucleon scattering length. A 
more accurate measurement of a+ not involving pion deuteron scattering data, as by the 



present PSI experiment R-98.01 [jrT|, will substantially reduce the uncertainty. 

We can therefore answer the question posed in the Introduction as follows: Chiral sym- 
metry in combination with already known parameters does not suffice to determine the pion 
deuteron scattering length due to two new contact terms with unknown couplings. However, 
once these terms are fixed from the pion deuteron scattering length they may be used to 
predict pion deuteron scattering data at non-zero momentum transfer. 



Acknowledgements 

We are grateful to S. Beane, T. Ericson, T. Hemmert, N. Kaiser, D. Phillips, A. Thomas, 
U. van Kolck and W. Weise for discussions and useful remarks. D. Gotta drew our attention 
to the ongoing experiments on pionic hydrogen at PSI. The INT (Seattle) provided generous 
hospitahty during part of this work (H.W.G.). We also acknowledge financial support by the 
DFG Sachbeihilfen GR 1887/1-1 and GR 1887/2-1 (H.W.G.) and by the Bundesministerium 
fiir Bildung und Forschung and the Deutsche Forschungsgemeinschaft. 



16 



A Calculating the One Body Scattering Amplitude 



As the pion can scatter off the proton or neutron inside the deuteron, the one body amphtude 
is (up to relativistic corrections) expected to be given by twice the physical amphtude A"^ 
only, irrespective of the presence of nested loops at NLO. We can confirm this formally 
and turn as an example to Fig. [l| (61). After the integration over the loop energy of the 
larger loop, the two loop integral to be performed is using the loop momentum assignment 
indicated in Fig. |l] (61) 



d^g d^ 



(27r)3 (27r)4 (g2 + ^ 



2^2 



Mir 



^2 _ ill _ /2 _ ^2)5 



(A.l) 



where the second term comes from the nucleon propagator inside the nested loop, and the 
last term is the pion propagator. Following the threshold expansion technique ||2^, we 
identify two poles in the energy integration of the nested loop: (1) M/q ~ g^, P ~ from 
the nucleon propagator, and (2) /q ~ |^| ~ '^n ~ Q from the pion. In the first case, the 
pion propagator is expanded because /q ^ 1^1? ^tt? so that the pion becomes instantaneous. 
Since the nucleon propagates forward in time, this contribution is expected to vanish. The 
dimensionally regularised loop integral over /q is indeed zero as no external scale is present 
in the denominators. In the second case, the pion pole is picked, and the nucleon propagator 
is expanded into 
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(A.2) 



The nucleon propagator becomes static as in HBxPT, and the two integrals factorise at 
LO in the expansion. The (relativistic) corrections are suppressed by powers of -p: ~ |, 
i.e. negligible in a NLO calculation. 

We therefore find (|3.3|) as conjectured to NLO, independent of any deuteron observable: 
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B Renormalising the Two Body Scattering Amplitude 

We will now construct the renormalisation group invariant (and hence physical) combina- 
tions of amplitudes in the two body sector of ird scattering at threshold. As mentioned 



(0) 



above, the parameter D2 enters in A^A^ scattering at NLO together with the parameter Cq 
which respects chiral symmetry. There, it is also needed to absorb the four dimensional UV 



divergence stemming from the one pion exchange diagram. We recall [|T8l that the pertinent 
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Figure 4: The contributions to NN scattering in the ^Si channel at NLO in ENT(KSW) 
which completely determine the renormalisation of 



terms to renormalise D2 in the triplet S wave channel of A^A^ scattering at NLO, depicted 
in Fig. ^, are in the centre-of-mass frame between nucleons with relative momentum p from 
one pion exchange 

■MA.^{p)y gl 
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and from the contact interactions involving D2 and C, 
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(B.l) 



^^.Tbte = - ( mlD2, bare(/i) + C^'LM 
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where Cq is the LO part of the strength of the two nucleon contact interaction Co, and 
the connected piece of the LO A^A^ scattering amplitude (lower line of Fig. ^) is 
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The one pion exchange part involving powers of the arbitrary regularisation parameter 
fi comes from the contact piece of one pion exchange which generates divergences in 3 
dimensions manifest in the PDS scheme. It is independent of the pion mass, i.e. unchanged 
by the chiral limit, and its regulator dependence can hence be absorbed into the definition 
of the renormalised coupling strength C^^^i •= ^o°bare + Ip accordance with the power 
counting On the other hand, the four dimensional divergence F giving rise to the 

logarithmic dependence on fi is chiral symmetry breaking and hence needs to be balanced 
by the definition of the bare two nucleon coupling -D2, bare- D2 is renormalised by setting 

2 
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At LO, Cq = — 37 (yU — 7) ^ is determined by demanding the triplet S wave to exhibit a 



pole at the deuteron binding energy in ( p.3|) . Therefore, one obtains finally 

^2,barc(/i)(/i - 7)^ 



9i 



l^r - 2 In 


'A* " 
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The combination of the divergent one pion exchange diagram and of the contact inter- 
action depending on the renormalisation group variant parameter D2, bare is traded for 
one, renormalisation group invariant parameter AJ^^. This dimension-ful, physical number 
parametrises the renormalisation group invariant strength of the chiral symmetry breaking 
contact interaction between nucleons which does not contain derivatives. It needs to be 
determined from experiment or from a microscopic calculation of A^A^ scattering in QCD. 
As only a variation of the pion mass can dis-entangle the effects of Cq°^ and D2 in A^A^ 
scattering at NLO, this process cannot serve to determine A^^^ experimentally. 

Inserting the definition ( [B.5| ) into the two body amplitude J?^^"^^ , (PT7|/|3.8|) reveals that 
some dependence on and on both logarithmic and power law divergences remains: D2 
serves only as a partial counter term to the diagrams Fig. (a) depending on and does 
not affect the divergence of the g\ diagram Fig. ^ ih). In the second, analogous step, the 
remaining divergences in the two body sector are easily shown to disappear by setting 
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where A^^ is another dimension-ful quantity, parameterising the renormahsation group in- 
variant strength of the quark mass independent couphng of two nucleons and two pions. 
The two body sector of the renormahsed amphtude of ird scattering at threshold, 
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6 V2 + 

contains therefore two un-determined parameters. Pion deuteron scattering at non-zero 
momentum transfer will allow to separate the two, and hence also Cq^^ from D2 in A^A^ 
scattering, but at present we choose to represent the two by one common scale 
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The final answer for the two body scattering amplitude in terms of renormalised quantities 
is therefore given by ( p.ll|) 
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The well known fact that dimensional regularisation (together with the employed renor- 
mahsation scheme) preserves chiral symmetry manifests itself in the observation that all 
appearing divergences are renormalised manifestly chirally invariant. 
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